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It is well established that the phase transition between survival and extinction in spreading models 
with short-range interactions is generically associated with the directed percolation (DP) universality 
class. In many realistic spreading processes, however, interactions are long ranged and well described 
by Levy-flights, i.e., by a probability distribution that decays in d dimensions with distance r as 
^-d-CT^ We employ the powerful methods of renormalized field theory to study DP with such long 
range, Levy-flight spreading in some depth. Our results unambiguously corroborate earlier findings 
that there are four renormalization group fixed points corresponding to, respectively, short-range 
Gaussian, Levy Gaussian, short-range DP and Levy DP, and that there are four lines in the {a, d) 
plane which separate the stability regions of these fixed points. When the stability line between 
short-range DP and Levy DP is crossed, all critical exponents change continuously. We calculate the 
exponents describing Levy DP to second order in e-expansion, and we compare our analytical results 
to the results of existing numerical simulations. Furthermore, we calculate the leading logarithmic 
corrections for several dynamical observables. 

PACS numbers: 64.60.ae, 64.60.ah,05.40.-a,64.60.Ht 
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I. INTRODUCTION 

The formation and the properties of random structures 
have been an exciting topic in statistical physics for many 
years. In the case that the formation of such structures 
obeys local rules, these processes can often be expressed 
in the language of epidemic spreading. It is well-known 
that two special spreading processes referred to in this 
language respectively as simple epidemic with recovery 
or Gribov process [l|, Q and epidemic with removal (gen- 
eral epidemic process) lead to random structures with 
the properties of percolation clusters: directed percola- 
tion 0, 0, @ in the former case and isotropic percolation 
in the latter. 

The Gribov process, also known in elementary parti- 
cle physics as Reggeon field theory (RFT) [1, 0? ], is 
a stochastic multiparticle process that describes the es- 
sential features of a vast number of growth phenomena of 
populations without exploitation of the environment near 
their extinction threshold. The transition between sur- 
vival and extinction of the population (infected individu- 
als) is a nonequilibrium continuous phase transition phe- 
nomenon and is characterized by universal scaling laws. 
The Gribov process belongs to the universality class of 
local growth processes with absor bing states [1, [13, [HI 
such as the contact process [l^ [l3l . Il4l | and certain cellu- 
lar automata [l^, , and it is relevant to a vast range of 
models in physics, chemistry, biology, and sociology. As 
usual, we refer to this universality class as the directed 
percolation (DP) universality class. For recent reviews 
see [11 [H. 

A continuum description of DP in terms of a den- 
sity n (r, t) of infected individuals typically arises from 
a coarse-graining procedure in which a large number of 



microscopic degrees of freedom are averaged out. Their 
influence is simply modelled as a Gaussian noise-term 
in a Langevin equation. The process has to respect the 
absorbing state condition: n (r, <) = is always a station- 
ary state. Then the minimal stochastic reaction-diffusion 
equation for the density n (r, t) is constructed as Q 



X~^dtn{r,t) ^ W^n{r,t) 
+ C(r,t) . 



T + —n{r,t) n (r, t) 



(1.1) 



The Gaussian noise ^ (r, t) must also respect the absorb- 
ing state condition, whence 



C (r, t) C (r', V) = X-\j' n (r, t) 5 {r - r') S [t - t') (1.2) 

up to subleading contributions. The history of the pro- 
cess in space and time defines directed percolation clus- 
ters in a ((i+ l)-dimcnsional space. The minimal process 
defined by Eqs. (|1.1|) and (|1.2p contains all the relevant 
terms needed for a proper field theoretic description of 
DP. 

In realistic situations the infection can be also long- 
ranged. One may think, e.g., of an orchard where fly- 
ing parasites contaminate the trees practically instanta- 
neous in a widespread manner if the timescale of the 
flights of the parasites is much shorter than the meso- 
scopic timescale of the epidemic process itself. Thus, 
following a suggestion of MoUison Q, Grassberger [g^l 
introduced a variation of the epidemic processes with an 
infection probability distribution P(r) which decays with 
the distance r as a power law, P(r) - r-'^-'^. We will 
somewhat casually refer to such long-range infection as 
Levy-flights although a true Levy-flight is deflned via its 
Fourier transform as P(q) ~ exp(— feg"^) with < cr < 2 
(to ensure positiveness of the distribution). 
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In Fourier space and in a long- wavelength expansion, 
the Langevin equation (jl.ip can be generalized to account 
for Levy- flights by a term proportional to q"n{q,t). In 
the case of 2 — cr = 2a > 0, the long- wavelength behavior 
is naively dominated by this new term. Grassberger cal- 
culated critical exponents in a 1-loop calculation which 
were discontinuous in the limit a +0, and therefore 
the applicability of the results was doubtful. In a for- 
mer paper [2l[, we have shown by applying the Wilson 
momentum shell renormalization group that only two of 
the critical exponents are independent in long-range DP, 
and that the critical exponents change continuously when 
the transition line between long-ranged and short-ranged 
spreading (with an a = ac < 0) is crossed. We have also 
shown that 



2P 

V 



1] 



(1.3) 



exactly, where /?, rj, and z are the usual exponents 
of short-ranged DP in d (transversal) dimensions. These 
results have been confirmed numerically by Hinrichsen 
and Howard [l^. Note that 77^ = z — 2 ~ rj, which cor- 
responds to the Fisher-exponent of equilibrium critical 
phenomena, is negativ here [23j . For a recent review for 
DP with long-range interactions see [23|. In this paper 
we reconsider the problem using methods of renormalized 
field theory in conjunction with an expansion in e and a. 

The remainder of this paper is organized as follows: 
Section |lT] reviews the field theoretic formulation of DP 
with Levy-flight spreading to set the stage, to provide 
background information and to establish notation. Sec- 
tion [n] first reviews the short-range limit of this model 
and then discusses our field theoretic analysis of the long- 
range limit. Section llVI represents the main part of this 
paper. It treats in detail the hybrid model for a = 0{e). 
Section |V] builds up on the results of Sec. llVl and present 
results for the critical exponents and logarithmic correc- 
tions of various dynamical observables. Section IVII con- 
tains some concluding remarks. 



II. MODELLING DP WITH LEVY-FLIGHT 
SPREADING 

To generalize the diffusional infection rate in the 
Langevin equation (|l.ip . we model spreading by writing 



9tn(r,t)|j„f = J d\'P{r-r')n{r\t) . (2.1) 

As a particular model for the positive distribution 
P(r) which contains all relevant properties, we use 
P(r) = PLR{r) + i-'sfl(r) with a short-range contri- 
bution Psfl(r) oc exp(— r^/a^) and a Levy-fiight part 
PLi?,(r) oc (r^ -I- 0^)"^''"'"°'"^. a denotes a microscopic 
length scale which, for simplicity, is assumed to be equal 
in both otherwise independent distributions. Fourier 
transformation leads to PsR{q) oc exp (— (a(7)^/4) and 



Plr{(1) oc K„/2{aq), where is the modified Bessel- 

function with index a/2. Long- wavelength expansion 
leads to 



P(0) - P(q) = A{aq)^ 



B 



[{aqY - {aqf 



+ 0(<^^g2+-) 



(2.2) 



with positive, non-singular, cr-dependent constants A 
and B. P(q) shows the two typical terms: diffusion 
^ [aq)^ and Levy-flights ^ {aqY ■ It is IR-stable, that is 
P(0) — P(q) is positive in both regions cr < 2 and a > 2 
if (7 ^ 0. Note the characteristic pole at cr = 2 that leads 
to a logarithmic contribution ^ ~{aqY \n{aq) to P{q) for 
CT — > 2. In the following we use for P{0) — P(q) its long- 
wavelength approximation (|2.2p . Whereas P{0) — -P(q) 
is always a positive real quantity, its long-wavelength ap- 
proximation changes the sign at a magnitude qg of the 
momentum of order aqg = 0(1) ifa > 2 or if cr < 2 
and A < B/{2 — a) which leads to a pole in the Greens 
function of the equation of motion. Of course this pole is 
a non-physical ghost; it arises when the long- wavelength 
approximation is used in a momentum regime where it is 
inapplicable. This happens in particular in dimensional 
regularization where integrations over internal momenta 
are extended to infinity. Hence, this method may be- 
come inconsistent if both types of g-dependencies, (ag)^ 
and {aqY , are used in common. We come back to this 
question in Sec. IIV CI In contrast, Wilsons momen- 
tum shell renormalization procedure avoids this danger- 
ous UV-region q — 0{l/a) because all momenta are re- 
stricted to the sphere q < 0(/i) with /i <C 1/a. This was 
the reason for using Wilson's renormalization group in 
our former publication [2lj. Note also, that dimensional 
regularization leads in the present problem like in other 
similar problems to the so-called triviality problem at the 
upper critical dimension, i.e., the dimensionally regular- 
ized theory misses logarithmic corrections and thus has 
to be viewed as an effective theory for low momenta. 

The stochastic equation of motion of the DP-process 
with the Levy-flight spreading and short-range diffusion 
can be written as 



X-'dtn{r,t) 



T+-gn{r,t) 



n (r, t) 



n{v,t)+C{r,t) , (2.3) 



where we set cr = 2(1 — a). Here the Levy- 
term on the right side is deflned in Fourier space as 
(— V^)^ "n(r,<) = g^^^^^^Ti (q, t) exp (iq • r). In or- 
der to develop a renormalized fleld theory, it is useful to 
recast the Lan gev in equation (j2.3p as a dynamic response 
functional |25|. [2611 



J [s, s] 



d'^rdtXSl A-^a 



+ ^(s-S) 



(2.4) 
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where s (r, t) ^ n (r, t) is the rescaled density which en- 
sures that g' = g and for which the time inversion symme- 
try s (r, t) ^ ~s (r, — t) (rapidity reversal in RFT) holds. 
s (r, t) is a response field that describes the response 
when a local particle source h (r, t) > is added to the 
Langevin equation ()2.3p . At the level of the dynamic re- 
sponse functional, this source leads to an additional term 
J d'^r dt h {r,t) s {r,t) in Eq. (|2.4p . Having the dynamic 
response functional, correlation and response functions 
can be computed as functional averages (path integrals) 
of monomials of s and s with weight exp {— Through- 
out this paper, functional integrals are interpreted in the 
sense of the so-called prepoint-discretization that sets the 
step function 9 (t) equal to zero for t = [131 • We stress 
that the usual short-ranged DP-model is recovered from 
the general expression of simply by setting c = 0, or 
c — 1 with a = 0. 

As a first step towards the renormalization group (RG) 
analysis of this model, we discuss its canonical scaling be- 
havior. Introducing the usual inverse length scale fi, we 
readily find s ^ s ^ /i''/^. For a > 0, the long-range 
Levy-term ^ (— V^) " naively dominates the usual dif- 
fusion term ~ V^. Hence, we may neglect the latter for 
a > 0, and we redefine (by rescaling of some parameters) 
c = 1. This produces an inverse time scale A/i'^, and 
T ~ /x'^ for the scaling of the control parameter. More- 
over, we obtain g^ where e = 2a — d = e — Aa (we 
will reserve the symbol e for the short-range case, i.e., 
e = 4 — d). The naive dimension of the coupling con- 
stant g allows us to identify the upper critical dimension 
dc{a) = 4 (1 — a) = 2a. This boundary separates trivial 
(mean-field or Gaussian) from non-trivial long-range be- 
havior if a > 0. Of course, the boundary a = 0, d > 4 
separates the regions with trivial long-range and trivial 
short-range DP, and the boundary a < 0, d ~ A separates 
trivial and non-trivial short-range DP. 



III. SHORT-RANGE AND LONG-RANGE 
MODELS 

We now turn to perturbation theory. In this section, 
we will first briefly review the short-range model obtained 
for c = 0, which has been discussed previously at many 
places (see [l^ and the references cited therein). Then 
we will treat, also briefly, the long-range model obtained 
for c — > oo. As usual in dynamical fleld theory, we focus 
on those correlation and response functions 



G 



N N 



(3.1) 



that require renormalization due to the presence of ul- 
traviolet (UV) divergences in Feynman diagrams as well 
as the corresponding one-particle irreducible (IPI) ver- 
tex functions with TV (A'^) external s- (s-) legs, F^^. For 
background on the methods of renormalization theory, 
we refer to 1281. 



A. The short-range model: c — 

We first review ordinary DP which is modeled by 
as given in Eq. (|2.4p with c = |9|, 0] . The upper criti- 
cal dimension is dc{0) = 4. Straightforward dimensional 
analysis shows that there are three superficially diver- 
gent vertex functions: Fi^i, ri^2 = —1^2,1, where the last 
relation follows from time inversion symmetry. In the fol- 
lowing, we use the superscript ° denote bare (unrenormal- 
ized) couplings, and we use the following renormalization 
scheme to cure the model of its UV divergences 

s^Z^/^s, l = Z^^^S, X = Z-^ZxX, 

f = Zl^ZrT + fc , g^ = G-^Z-^Zl^ZuU[i'' , (3.2) 

where = F(l + £/2)/(47r)''/2 

is a convenient ampli- 
tude, u represents the dimensionless coupling constant, 
and the control parameter r is zero at the critical point. 
In dimensional regularization, the critical bare value of 
the control parameter, Tc, is of the form 



°g^''S[e), 



(3.3) 



where the Symanzik function S(e) has simple IR-poles 
at each e = 2/fc with fc = 1, 2, .... Hence, fc is not 
a perturbational quantity and is formally zero in the e- 
expansion. Note, however, that minimal renormalization, 
i.e., dimensional regularization in conjunction with min- 
imal subtraction, does not imply the e-expansion [2^ . 
The renormalization factors Z... are functions of u and 
have in minimal renormalization the expansions 



y(")(M) 



Y 



-u . 



l—n 



I 



(3.4) 

where the Z... are determined in such a way that the 
perturbation expansions of renormalized quantities are 
free of singularities if e goes to zero. They are given to 
second order by [§, [l^, [l^ 

4e \e 2 3 / 32£ ^ ' ' 



u 

8£ 

u 
2e 
2u 



31 


35 


" T 


+ y 














I) 





u 



3 J 128i 



+ O {u^) , 



(3.5) 



A renormalization group equation (RGE) for the model 
can be derived in a routine fashion by exploiting the fact 
that the unrenormalized response and correlation func- 
tions have to be independent of the inverse length scale 
/i introduced by renormalization. This reasoning leads 
straightforwardly to the RGE 



_ N + N 
— 7, — 7 



^N,N ~ ^ ' 



(3.6) 
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with an RGE differential operator V = l^d/dfi\hiiio (the 
I bare indicates that bare quantities are kept fixed while 
taking the derivates) given by 

The RG functions result from the finite logarithmic 
derivatives of the renormalization factors 



7... = 



d\nZ 



d 



bare 



1=1 



(3.8) 



/3k = M 

C= 

K = fl 



du 



a In A 



bare 



51nT 



(-£ + 27a + 7-7u)u, 



(3.9) 



Their perturbation expansions are 



0«^-^+(2-31„l)5^ + O(„3) 



<(«) 



u 

8 



4\ li2 

17 -2 In- + 0(u^ 

3 / 256 ^ 



.(.)^-- 7+lOln- - + 0(.3 



/34?i) 



3u / 4\ 
-^ + Y-^«^ + 10«1^3 128 



(3.10) 



The asymptotic solution of the RGE. Eq. (|3.6p . leads to 
the stable fixed point u = with given by 



£)p, , 2£ 

= (e) = — 



1+ ( l69 + lO61n0 ^ + 0(e2) 



(3.11) 

as the stable solution of f3u{ut,) — 0, and to the scaling 
form 

G^,^({r,t},r) = /(^+*H'^+"-'/2 

xG^,v({'r,r-«0,ri/''^-T) (3.12) 

of the response and correlation functions, with the three 
independent critical exponents 

VSR = l{u^^.) , l/ysR = '^- i^{u^) , zsR^2 + ({u^) . 

(3.13) 

These are the very well known critical exponent for short- 
range DP. To second order in e-cxpansion, they are given 



by 



VSR 



zsR^2-- 



25 161, 4 , ^, 2n 
288 + 144 ^^3^^ + ^(^) 



67 59 , 4\ ^. 2, 

\ In- e + 0(s^^ 

288 144 3 ' ^ 



l^SR 



16 



1+r^ 



17 4^ 

— In - 1 e + 0(e' 
V288 144 3.' ^ 



(3.14) 



B. The long-range model: c — » oo 

As we have shown in [2l| by using Wilson's 
momentum-shell renormalization group, and as we dis- 
cussed in the introduction, the discontinuity of short- 
range and long-range critical exponents at a = is spu- 
rious and can be remedied. In renormalized field the- 
ory, the key is to recognize [30| that there is a region 
of small a = 0(e), where a careful analysis of the RG 
flow reveals a smooth connection between the a < 0, 
a = 0(e), and a > regions. Here, we analyze the last 
case, which belongs to the true long-range region. The 
case a = 0(e), being a "hybrid" between short-range and 
long-range models, will be deferred to the next subsec- 
tion. 

We recall from our discussion at the end of Sec. |TT] 
that the upper critical dimension for a > is dc{a) — 
4 (1 — a) = 2(7, and we define e = 2a — d, to be distin- 
guished from e = 4 — (i. Considering the response func- 
tional J, Eq. (12. 4p with c 7^ 0, we see that the operator 

sV^s is superficially irrelevant compared to s (— V^) ^ " s 
and may be dropped formally in the limit c 00. This 
limit is feasible after the rescaling A — > A/c, r — *■ re, and 
g gc. The canonical dimensions of the fields do not 
change compared to the short range model. However, 
the inverse time-scale changes to Ay^t*^, and the canonical 
dimensions of the remaining parameters are t ^ fi"^ and 
g ~ /x'^''^. As above Fi.iand ri.2 = — r2,i are superfi- 
cially divergent for e — s- 0. Moreover, all divergent con- 
tributions to any vertex function are polynomial in the 
momenta, so that the operator s V^) " s needs no 
counterterm. Hence, we use the renormalization scheme 

T = ZrT + fc, g^ = Ag^ Z-^ ZuU^^ , (3.15) 
which produces the renormalized response functional 



d'^r dt As 



X-'Zdt 



ZrT+i-^^y 



(3.16) 



Here and in the following we use an overbar to distinguish 
the renormalization factors of the long-range and hybrid 
models from those of the short range model. Ag is a 
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suitable amplitude whose precise definition will be given 
later. Here, in minimal renormalization, = g^"^^^ S{e) 
with an appropriate Symanzik function 5(e) having sim- 
ple poles at e = cr/fc with fc = 1, 2, .... Note by com- 
paring Eq. (|3.2p with Eq. (|3.15p that the renormalization 
schemes for the short-range and the long-range models 
are of the same form except for the renormalization fac- 
tor of the kinetic coefficient, which now is Z\ — 1. Here, 
we have the expansions 



Z = 1 



n=l 



Y 



(") 



I 

' u 



n=l 



(3.17) 

Of course, the functions y(")(u) are different from the 
functions Y^^\u) and have be determined by perturba- 
tion theory. Nevertheless, the RG functions for the long- 
range case can be transcribed from the short-range case, 
Eqs. dS^, ([XT]) and simply by decorating each RG- 
hmction with an ovcrbar and setting 7a = 0, i.e.. 



- du 



bare 

9 In A 



9/i 



= (-£ + 7-7«)m, 

91nT 



bare 



(3.18) 



The scaling form of the response and correlation func- 
tions follows from the RGE as 

G^,^({r,0,r)=/(^+*)('^+''-)/2 

xG^^^({Zr,P-«t},/-i/''-«r), (3.19) 

with the rjLii and v^ji given by 

= 7("*) , l/i^Lfi = cr - k{u^) . (3.20) 

The third exponent, zlr, is related to these exponents 
by the exact relation 



ZLR 



VLR ■ 



(3.21) 



Hence, we have only two independent critical exponents 
in the long-range case, viz. tilr and vlr- 

After this general discussion, we now turn to our actual 
perturbation calculation. The propagator of the theory 
reads G{q,t) = 9{t) exp {—X{t + \q\'^)t) in wavevector- 
time representation. Note that the prepoint discretiza- 
tion mandates that we have to use 9{0) = through- 
out [131. Note also, that the propagator is free of the 
ghost-problem, i.e., the exponent is always negative. The 
1-loop self-energy with frequency cu and wavevector q as 
the first contribution to the vertex function Fi i reads 



1 



it^/A + 2r+ |p + q/2r 



q/2r ■ 

(3.22) 



To calculate this self-energy, it is useful to expand in q 
and Lo and to use the identity 



2^-'^/2r(i + d/2) J d^pf {\pn 

= 27r-'^/-r(l + d/a) J d''/yf{\p\^) 



(3.23) 



which leads the calculation of the primitive divergent 1- 
loop diagrams back to integrals of the usual known short- 
range type. We obtain 



At > + finite , 



^q,.) = -f;.-/^^.{..4-^A.} 

(3.24) 

where we have displayed only the pole-terms in e. The 
pole at e = cr is an IR-pole and can be removed by in- 
troducing a new mass parameter m fi instead of t via 
introducing m""' = d\nri^i{Lj ^ 0,q,T)/dq'^\^^^ (keep- 
ing in mind that it is an IR-poles, one may also simply 
ignore it in the e-expansion) . The amplitude Ag is de- 
fined by 



A, = 



r (2 ^ e/g) r (1 -f e/a) 
r(cr- e/2) (47r) 



d/2 



(3.25) 



Note, that Ag becomes if tr 2. Expanding 
{^'^ / tY / / {\ — e/tj) in e, and using the renormalization 
scheme Eq. p.isp . we arrive at the singular part of the 
vertex function Fi i in 1-loop approximation 

ri,i (q, uj) ^ {ZtLj + ZrXr) - - I — + — 

(3.26) 

As announced above, there is no singular term propor- 
tional to |q|'^. Using the same techniques, wc get 



Ti,: 



?l/2 



(3.27) 



for the singular part of the other superficially divergent 
vertex functions. We read off the 1-loop renormalizations 
of the long-range model 

Z=l + ^ + 0{u^) , Z, = l + ^ + 0{u^) , 

Zu = l + ^ + 0{u^) . (3.28) 

The RG functions here have the expansions 

C = 7 = -J" + 0(m2), S=iu + 0(u2), 



-e+ ^u + 0(w^) ) u, 



(3.29) 



instead of those given in Eq. (|3.18p for the short-range 
model. The stable fixed-point is m, = u^^ = 4e/7 -I- 
O(e^), and the expansions to first order of the long-range 
critical exponents are 

VLR = + (e^) , i/^^j, = a-y + O (e^) , (3.30) 
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which should be compared to the short-range exponents 
given in Eq. (|3.20p . As has to be the case, our 1-loop 
results p.30p are in perfect agreement with the 1-loop 
results for the long-range exponents derived in [2]| by 
Wilson's method. 



IV. THE HYBRID MODEL: a = 0(e) 

Here, we turn to the analysis of the key region in (d, a) 
space, namely, a = 0(e). The naive a — > limit of 
the long-range model presupposes e a and hence fails 
to resolve the crossover between the SR and LR models 
which occurs for a = 0{e). For both a and e small, we 
follow the work of Honkonen and Nalimov [301 ■ 



of the renormalization parts of the diagrams are tamed 
by lower order counterterms. After time integrations over 
the {V — 1) time segments each between vertices, any di- 
agram has the qualitative form 



— 9 



k=0 
oo 



{m? -t- p^) 



[d'^p)^ 



OO 

E 

fc=0 



n \A(L.V,k) k-V 

CL,v,kA ' 'c g . 



(4.4) 



A. Renormalization 

Our starting point is the renormalized response func- 
tional 

+ [z.f-Z,V2 + c(-V2) (4.1) 
where we use the renormalization scheme 



\ ~ Z Z\X , 



f = Z-'Zrf + g'^ G-'Z-'Z^'Z^ufi' , 

c - Z^'i V" , (4.2) 



and the abbreviations g = \/uii^^^ and c = wfi" 

" s does not need a counter 



.i^". As 

before, the term s (— V^) 
term as long as a 7^ 0. 

Using the approach by he Honkonen and Nalimov, we 
construct our renormalization factors Z... by generaliz- 
ing the minimal renormalzation program that led to the 
expansions p.4p and (|3.17p . Here, the renormalization 
factors are now functions of u and w, and they contain 
poles of all linear combinations Si^k = le + 2ka with / — 1, 
2, . . . and fc = 0, 1, 2, . . .: 



Z 



00 oc 

EE 

1=1 k=0 



Y 



(1) 



■.Lk 



le + 2ka 



(4.3) 



where the coefficients F*"^; f. can be chosen such that they 
are independent of s and a if both are of the same order. 
To explain this, let us consider the primitive divergence 
of an irreducible diagram consisting of V vertices P prop- 
agators L independent loops, and E external amputated 
legs. By definition, a primitive divergence is a divergence 
that arises if all the Ld inner momentum integrations 
tend uniformly to infinity. In primitive diagrams, they 
are the only UV-divergences. In non-primitive diagrams, 
they are the divergences that remain after all divergences 



Here, the mass is a linear combination of the con- 
trol parameter r and frequencies uj and serves as an IR 
regulator. A is a momentum cutoff, and A{L,V,k) ~ 
dL — 2{V — 1) — 2ka is the degree of primitive diver- 
gence of the diagram. Using the topological relations 
P = L+{V-1) and 3y = 2P+E, hence V = {E-2)+2L, 
we obtain 



A(L, U, k) = 2(3 -E)- {Le + 2ka) 



(4.5) 



The first part, 2(3 — E'), is the superficial divergence of 
the diagram with E legs. The second part, — [Le + 2ka), 
denotes the combination which is converted to a simple 
pole in dimensional regularization. This pole-term must 
be eliminated by a counterterm if both e and a become 
small quantities. Thus, the overall form of the renormal- 
ization factor is as given in Eq. (|4.3[) . Up to now, the 

constants j, may still be functions of e and a. How- 
ever, if a/s is finite in the limit e — > 0, we can neglect 
this dependencies in the sense of the minimal renormal- 
ization. Hence, we can apply this Honkonen-Nalimov 
scheme only if a — 0{e). 

Next, we calculate the logarithmic derivatives of the 
renormalization factors. Note that in minimal renor- 
malization the only terms of the /3-functions, [3^ = 
/i9u/(9/i|baro and (3^ = fJ-dw/d^lhaici which contain e and 
a explicitly come from the //-factors making g^ and c di- 
mensionless, (cf. Eq. (|4.2[) ). i.e., (3u = —eu + ... and 
(3w = ^2aw + • • ■ . Thus, we obtain from Eq. 



7... = M- 



dlnZ 



dfj. 



d 



^ ^ 00 00 Y 



' ^ 9u ^ dw ) ^ ^ le -I- 2ka 

l=\ k=0 



(1) 



EE- 



■ ■■■,l,k k-l 
W U 



00 00 



EE^^ 

1=1 k=0 



(4.6) 



which should be compared to Eq. (|3.8p of the short-range 
case. 
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Now wc can the relate the functions 7...(?2, w) to the 
original DP-functions 7...('^), Eq- P-8p . pertaining to the 
short-range case. To this end, we consider the hybrid 
response functional and its renormalizations for a = 0: 

^[s,s]U^o== / dtd'^rXslzX-^dt + Zrf- {Zx + w)V^ 



(4.7) 



This functional takes on the form of the original short- 
range DP response functional, 

J[s,s]= J dtd'^r Xsl^ZX-^dt + ZrT ~ Zx\7^ 



1/29 



(s-s) }s, 



(4- 



if we identify the parameters 



X={l + w)X, T = (l + w) T, u = {1 + w) u 

(4.9) 

and the renormalization factors 

Z (u) ~ Z (w, w) , {1 + w) Zx {u) = Zx (m, w) + w , 

Zr (u) = Zr {u, w) , Zu (u) = Z.^ {u, w) . (4.10) 

The last identifications lead, by comparison of Eq. (|4.3p 
with Eq. (|3.4p . to the relations 



(4.11) 



in the case of Z, Z^, and Zu, and to 



E^i;iV- = r-/(i + ^) 

k=0 



1-21 



(4.12) 



in the case of Zx- Collecting, we obtain for the logarith- 
mic derivatives, Eq. 



■^(m,^) = 7(u) , 

Jr {U, W) = 7r (u) , 

7„ (w, w) = 7„ (u) , 

7A iu,w) = {1 + w) 7A (u) 



(4.13) 



Using the renormalizations (|4.2p . the RG functions be- 
come, 



- ait 

Pa = 



= (-£ + 27a + 7-7„)m. 



dw 



= (-2a + 7A) w. 



bare 

a In A 



dfj. 
91nf 



bare 



K = /i- 



= 7A - 7r 



(4.14) 



However, Eqs. (|4.9p and (|4.13p suggest that it is more ap- 
propriate to use instead of A, f, u, and w the parameters 
A, T, and u defined by Eq. (|4.9p . and 



2aw 

1 + w 



(4.15) 



as the parameters of the theory. The response functional 
then becomes 

J [s, s] = J dtd'^rXsl^ZX-'^dt + ZrT - ZaV^ 



2a 



'y^f(^-s) 



(4.16) 



Note, that J' coincides up to a rcscaling with J7l_r, 
Eq. (|3.16p . if V ~ 2a. Note also, that the gradient terms 
of this response functional take the same form as pro- 
posed for the long-wavelength or gradient expansion of 
the general long-range spreading, Eq. (|2.2p . Hence, we 
expect the same difficulties concerning the positivity of 
the propagator for momenta q > fj. and cr > 2 if i; > 2q;. 

The RG functions for the new variables are easily de- 
rived from Eq. (j4.14p by using Eqs. (|4.13p . 



/5« = M 



du 
9/i 



= (-£ + 2v + 27a + 7 - 7«) w : 



bare 



K = ^ 



a In A 



dfj, 
ainr 



(—2a + V + 7a) V , 



= 1 - 1\ - V , 



1\- It + V . 



(4.17) 



bare 



In comparison with the short-range RG functions, 
Eq. (|3.9p . the new GcU-Mann-Low function (3y as well as 
the functions /3„, Q, and k have only an additive contri- 
bution of the variable v. Note that we have to consider 
both u and v as being parameters of order £ ~ a. In 
terms of the new variables, the RGE is still of the form 
given in Eq. p.6p . however, here 



9 .^ C) d 

o^i oX OT 



-0 — + /3 — 
ou ov 



(4.18) 



has to be inserted as the RG differential operator. 



B. Asymptotic scaling regions in (e, Q)-expansion 

In this subsection wc analyze the different scaling re- 
gions in a (d, cr)-diagram using the results on the hybrid 
model derived in Sec. IIV Al The asymptotic scaling of 
response and correlation functions is governed by the var- 
ious fixed points of the renormalization group. To find 
stable fixed points of the RGE and the corresponding 
scaling behavior, we have to find solutions of the equa- 
tions (3u = Pv = with the GcU-Mann-Low functions /?„ 
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and j3y as given in Eqs. (|4.17p . The different fixed points 
can be elassified by setting to zero the different factors 
of these functions. 

There is tlie trivial fixed point 0. To find 

its stabihty conditions we determine the eigenvalues of 
its stability matrix 



/3* 



(dl3u/dul d/3u/dv\ 

[dp^/dui dp,,/dvi 



— £ 







-2a 



(4.19) 



The eigenvalues of this matrix, o^i = — e and lu2 = —2a, 
are positive for d > 4 and cr > 2, i.e., we retrieve the 
mean-field region of short-range DP. 

Another fixed point, the trivial long-range fixed point, 
is given by m* = 0, = 2a. Its stability matrix reads 



/4a - e 
a - l,2a7L 2a 



(4.20) 



where the stroke at denotes the derivative with respect 
to u. The eigenvalues uji = 4a — e and lu2 ~ 2a arc 
positive for d > 4(1 — a) — 2a, a < 2 which marks the 
region of stability of the trivial long-range fixed point. 

Next we come to the fixed point = 0, > 0, as the 
solution of 27>,, -I- 7* — 7ti* = £• Of course, this is the 
fixed point of the normal short-range DP with 
The stability matrix is 



,DP 



7a* - 2a 



(4.21) 



The first eigenvalue of this matrix, uji^{2j'^^ + 7* ^ 
7u*)'"* = £ + 0(£^) shows the stability range of non- 
trivial short-range DP: d < 4. Using 7^, — 7a (w?^) — 
^SR _|_ 2 — z^^, wc find the stability condition against 
long-range spreading: 



a > z 



SR 



SR 



(4.22) 



Now, we come to the interesting LR-region where > 
and w* 7^ 0. In this domain, the stable fixed points of 
(|4.17p are solutions of the fixed point equations 



2a - 7A* > 



4a — 7* 



lu 



(4.23) 



Using the 7- functions which follow from Eq. (|3.5p . and 
which have been utilized in Eq. (j3.10[) . wc find the fixed 
point 



,LR. 



4£ 

y 



1 50 + 9 In ■ 



3/98 



O (£^a£) 



= 2a 



£ 

14 



1 - 



17- 526 In - 
3 



£ 

392 



0{e\ae] 



(4.24) 



The critical exponents in the LR-region are found from 
Eq. gUl) as TiLR = 7("* ^) and 1/ulr = 2- k{uI:^) 
with the expansions 



= [l + c„(a)£H-0(£2)] , 
I/VLR = - y [1 + c^(a)£ -I- O (e^)] , 



(4.25) 



where 



c^(a) 



36, 4 
— In - 

49 3 



9 , 4 

1 In - 

98 98 3 



0{a) , 
0(a) . 



(4.26) 



Note that wc have already encountered the first order 
contributions in Eq. (|3.30|) above. Of course, z^n follows 
from the exact relation p.2ip . The stability matrix for 
this fixed point reads 



R — ( 2m, 
= V7a,«* ^* 



(4.27) 



with A = (27;, +i:-lL) = 3/2 + 0{e) > 0. The two 
eigenvalues arc given by 



UJ± = 




- {A-2y^Ju,iu 



(4.28) 

They arc positive as long as (A — 27^^) m*w* = (7^ — 
7^^)m*w* = (49/16 + 0(£))£v* > 0. This condition leads 
to £ > and z;* = 2a — 7a* > 0. The long-range fixed 
point looses its stability if the line 2a = 7a* is reached. 
At that point, all critical exponents change over contin- 
uously to the usual short-range DP-exponents as can be 
easily seen from Eq. (|4.17p . Hence, stability boundary of 
the long-range Levy-flight exponent c = 2(1 — a) is given 
by 



ZSR - VSR 



12 



17 

288 



— In- )£ + 0(£2,a£) 



(4.29) 



which is less than two. This fact is astonishing because 
for cr lower then 2 but greater then a^, the long-range 
part q'^ is naively irrelevant in comparison to the normal 
diffusional part q^. However, in an interacting theory it 
is not the free propagator but rather the response func- 



tion x(w, q; r) = Pi, i(w, q; r) ^ = q"! 
which is the deciding quantity. Hence, 



f{u/q^,T/q^/-) 
one has to com- 



pare q" with q^ 



to find out which is leading for 



g — > 0. The other stability boundary is approached if 
u* goes to zero, and is given by £ = £ — 4a = 0. This 
boundary coincides with the value found above: d ~ 2a. 
At this line the exponents cross over to their long-range 
mean-field values. 

To summarize our findings regarding the scaling re- 
gions: all boundaries between the four scaling regions, 
namely short-ranged DP, long-ranged DP, as well their 
two mean- field counterparts, are generally given by the 
four lines in a [d, cr)-diagram where one or two of the 
fixed point values u* and t;* vanish. There is no room for 
other stability lines as some authors argued [3l| . 
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C. Landau's ghost 

As we have remarked at several points, the propagator 
becomes problematic for higher momenta. Consider the 
g-dependent part of the inverse rcnormalizcd propagator 
of the hybrid theory, 



G(q,0,0)-i = Ag2 



1 



2a 



-2a 



1 



(4.30) 



where u > to ensure stability (positivity) for q s- 
for both signs of a. The part proportional to v changes 
sign and leads to a loss of stability for negative a at a 
momentum qg given by 

Hqg/^^) = 2^ ln(l + 2\a\/v) ^ 1/v (4.31) 

for small a. This ghost reminds of Landau's ghost in 
quantum electrodynamics. The momentum of the ghost 
goes exponentially to infinity if w — > 0. This ghost even 
arises for positive a if w > 2a. The fixed point of the 
hybrid theory does always belong to this region! The 
correct interpretation is the following: our asymptotic 
theory is just an effective theory in the sense that it can 
only be used in the low momentum limit in a perturbation 
expansion with v as an expansion parameter. Therefore, 
the second part of Eq. (|4.30p must be considered as a 
perturbation. 

Let us demonstrate this in some detail for the inverse 
response function Fi^i as an example. To this end, we 
work to 1-loop order using u, v, e, and a as first order 
quantities. The zeroth order is 



ri°l(q,- 



0,r = 0) = Xq^ 



(4.32) 



Using renormalized perturbation theory, adding both 
first order terms, and neglecting higher order terms, we 
obtain 



As expected, the RG proves to be the systematic tool to 
resum all logarithms to yield the correct critical expo- 
nent. 

This procedure holds for all a = O (e) irrespective 
of the sign as long as the fixed point with > is 
stable. Otherwise ~ 0, — 2e/3, and one gets 
Fi 1 ~ Ag^^''-5« with risjf = — e/12, i.e., the known be- 
havior for the short-range case. We once more point out 
that the second term (the Levy-flight contribution) has 
to be handled as a perturbation to the desired order and 
not as a part of the unperturbed propagator. Also, one 
has to interpret the special case a = as a relevant log- 
arithmic perturbation ~ vq^ln^ji/q) of g^. Only if v is 
strictly zero the short range case is recovered. Thus, one 
cannot expect a continuous behavior at a = comparing 
the short-range versus the Levy-flight directed percola- 
tion. 



V. CRITICAL BEHAVIOR OF DYNAMIC 
OBSERVABLES IN LEVY-FLIGHT DIRECTED 
PERCOLATION 

In this section we will harvest some of our previous 
results to calculate scaling forms and logarithmic cor- 
rection for those dynamic quantities in long-ranged DP 
that are most suitable from the vantage point of nu- 
merical simulations [2^ . Two key observables with re- 
spect to simulations are the density of infected individu- 
als p{t) = {s{r,t))pg for t > if the initial state at time 
i = is prepared with a homogeneous initial density po, 
and the response function x(r,i) = (s(r, t)s(0, 0)) that 
yields the density of infected individuals after the epi- 
demic is initialized by a pointlike source at t — and 
r = 0. 



A. Scaling properties 



ri,i(q,o,o) 

= Xq' 




1 - - In 2 Aq 



^In A 
8 \2p 



(4.33) 



Now we use the fixed point result = 2a — ^x* with 
7A = -u/8, to get 



ri,i(q,c^ = 0,T = 0) 
^(l-!^ln2)A.^ 



1 - 2aln 



(4.34) 



The scaling properties of the density of infected indi- 
viduals and the response function follow from the RGE 
p.6p taken at the long-range fixed point of Eq. (|4.17p and 
by identifying x(r, t) and p{t) with the Green functions 
Gi.i(r, i;T) and G'i_o(r, t, po), respectively. The initial 
density po is introduced into the the response functional 
via a (bare) source h{r,t) = po5{t) with renormalization 
Po = Z~^^^po [1^, which leads to an additional deriva- 
tive term ^■ypod/dpo in the RGE. We obtain the scaling 
forms 



p(^)=^-^5p(r^l/^^po^'+'), 



(5.1) 



where the 5... are appropriate scaling functions. We drop 
in this section all the subscripts at the critical exponents 
because we are interested in the long-range case only. 
Hence, 77 = tilr and 1/z/ = 1/vlr, see Eqs. (|4.25|) . and 
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z = a + r], 5={d + ri)/2z, 9 = —rj/z. The expansions of 
the latter two are 

S = l-^[l + csia)e + 0{s^)] , 
e^^[l + cg{a)e + 0{e^)] , (5.2) 

where 

CA (a) = I In — I + O (a) , 

^ ' \294: 49 3 J ^ ' ' 

At the critical point t = 0, Eqs. (|5.ip show that the 
mean square radius of spreading from the origin scales as 
average number of infected indi- 
viduals N{t) ~ J d'^rx{r,t) ^ . Starting with a homo- 
geneous finite value po, the critical density first increases 
in a universal time regime with the same exponent as 
pit) ~ Pot^ ■ Then, after some crossover time, it decreases 
as p{t) ~ . If one starts with a full lattice of infected 
sites corresponding to an infinite initial value poi only 
the last scaling behavior is seen. Because of asymptotic 
time-reflection invariance of DP (duality symmetry) this 
behavior characterizes also the survival probability [33 |. 
P{t) - t-^. 

To first order all exponents arc identical to the LR- 
exponents derived from Eqs. p.30p . of course. The full 
Qf-content of the functions c...(q;) in Eqs. (|4.25[) and (|5.2[) 
must be calculated from the long-range model, Eq. (|3.16p . 
Hence, one has to be careful when applying the ex- 
pansions to O (e^) for a = 0, e.g., in d = 1 where 
a~ (3 — e)/4. In Figs.[l]to[3l we compare our results for 
the exponents z, (5 and 9 to numerical results for d — 1 
by Hinrichsen [l8|. For this comparison, we use the first 
order expansions in e (red curves), which are exact in a, 
and the second order expansions in e where we neglect the 
a dependent parts of the c...(a) (green and blue curves). 
The green curves show our second order result for z and 
results obtained for 8 and 9 by using without further ex- 
pansion the scaling relations relating 5 and 9 to z and 77. 
The blue curves stem from using these scaling relations 
and then properly expanding 5 and 9 to second order in e. 
For (T in the range from 1/2 to roughly 1, the numerical 
data and the analytic results agree remarkably well. For 
larger a sigma, the agreement suffers, but is well within 
the expectations for the methods used here. 

B. Logarithmic corrections 

Above the boundary between the genuine and the triv- 
ial long-range regions (d > 2cr, ct < 2), the couphng con- 
stant g tends to zero under the RG. However, g represents 
a dangerously irrelevant variable here, since it scales var- 
ious observables, and setting 5 = rigorously leads ci- 
ther to zero or infinity for relevant quantities. Due to 



Z 




0.5 1.0 1.5 2.0 

a 

FIG. 1: (Color online) The exponent z as a function of a for 
d = 1. The data points stem from simulations by data by 
Hinrichsen [l^- The red (upper) and the green (lower) curve 
correspond to our 1-loop and 2-loop results, respectively. 




0.5 1.0 1.5 2.0 

O" 

FIG. 2: (Color online) The exponent (5 as a function of a for 
d = 1. For an explanation of the curves, see the main text. 




0.5 1.0 1.5 2.0 

a 

FIG. 3: (Color online) The exponent S as a function of a for 
d = 1. The color-coding of the curves is the same as in Fig. [2] 
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its twofold nature as both a relevant scaling variable and 
an irrelevant loop-expansion generating parameter, g has 
to be treated very carefully. To set the stage for such a 
treatment, let us briefly review a few fundamentals of 
dynamic field theory. In broad terms, one attempts to 
determine the cumulant generating functional defined by 
the functional integral 



Wlr H,H =ln V [S, s] exp - Jlr [5, s] + (ff, s 



+ (H, S) 



(5.4) 



Functional derivatives with respect to the sources H and 
H define the Green's functions. The generating func- 
tional for the vertex functions Flr [s, s], the dynamic free 
energy, is related to the cumulant generating functional 
via the Lcgcndrc transformation 



Tlr [s, s] + Wlr 



(5.5) 



with s = SWi,B./SH and s ~ SWi,b./SH, and vice versa. 
In terms of Flr, the twofold nature of g is lucidly exposed 
by writing 

rLR[s, s; T, g] = g"^$LR[5s, gs; r, u] . (5.6) 

The expansion of the functional <I>lr[(7s, gs; t, u] into a 
series with respect to u yields the loop expansion. The 
zeroth term _g~^$LR[.9^i 9^'^ ''"i 0] is just the response func- 
tional JTlr, Eq. (|3.16p . itself. The scaling form of the gen- 
erating functional for the cumulants that corresponds to 
Eq. (1121) reads 

WMH,H;T,g]^g-^ni^R[gH,gH;T,u]. (5.7) 

To leading order in the logarithmic corrections, we may 
neglect the dependence of fl and $ on u. Functional 
derivation lead the to the Green's functions 

G^.^({r, t}, r; u) ^ u-'+^^+'^^/^T^^^ ({r, t}, r, u'/'po 

(5.8) 

where Tj^ ^ arc the contributions of loopless trees con- 
sisting of iV + TV — 1 propagators and iV -I- iV — 2 vertices. 

The characteristic equations that follow from the RG- 
functions, Eqs. (|3.29p . with p,{£) = Ifji and d = 2(7, are to 
lowest order given by 



M 


= PW)) = 




d\nX{e) 


= 7(i2(£)) = 


4 


d\nXx{(.) 
dliil 


= ~Q{m) = 


4 


d\nXr[(.) 
d\nt 


= R{u{^)) = 


2 ' 



(5.9) 



Solving these equations, we obtain asymptotically for 
1, 



|ln^P 
|ln£| 



1/7 



-2/7 



(5.10) 



X T, 



X T 



Hence, we get 

GAr,7v({l-,0,-r;w) 

^u{i)-^[u{i)e"x{t)\^''^'^^'^ 

({^r, rxA(€)t}, i-"Xr{i)T, u(€)i/2r"x(^)i/2po 

lln^l riln^r^/^ 

({^r, r |ln^|i/^ i}, r'^ |ln£r'/^ r, r'^ \\n£\-^'^ po) 

(5.11) 

as solutions of the entire RGE . Choosing either 
r |ln^|^/^i ~ 1 or |ln^|"^/^r ~ 1, we deduce that at 
the critical point 

p{t)^po{\ntf" (5.12) 
in the initial time region, 

p(t)-P(i)-<-i (Int)^/^ (5.13) 
in the late time region, and 

2/(j 



(5.14) 



VI. CONCLUDING REMARKS 



In summary, we studied DP with Levy-flight spreading 
by using the powerful methods of renormalized field the- 
ory. Our work confirms the previously known RG fixed 
point structure including their stability regions and the 
fact that the critical exponents change continuously in 
the crossover between short-range DP and Levy DP. We 
calculated the critical exponents for Levy DP, which have 
hitherto been known to first order, to second order in an 
expansion in e and a. These results agree well with the 
existing numerical simulations for d = 1. In addition, 
we calculated the leading logarithmic corrections for sev- 
eral dynamical obscrvables that are typically measured 
in simulations. 

We hope that our work stimulates further interest in 
long-range DP. It would be interesting to see further sim- 
ulation results, e.g., for the critical exponents for d > 1 
and for logarithmic corrections. Also, it would be in- 
teresting to have analytical and numerical results for 
other universal quantities such as scaling functions and 
amplitudes. In a forthcoming paper we will apply the 
same methods to the long-range GEP, that is to dynamic 
isotropic percolation. 
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